Berry Phase and Spin Qunatum Hall Effect in the Vortex State of Superfluid 3 He in 

two dimensions 
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We show that the spin quantum Hall effect in the vortex state of two-dimensional rotating su- 
perfluid 3 He can be described as an adiabatic spin transport of Bloch quasiparticles. We show that 
the spin Hall conductivity is written by the Berry phase as well as the Chern number. The results 
have similarity to the adiabatic pumping of Bloch electrons and the spontaneous polarization in 
crystalline dielectrics. 



I. INTRODUCTION 



The Berry phase (the geometrical phase) arises in 
quantum mechanical systems with an adiabatic change 
on a closed loop in a parameter space Q . In spite of the 
fact that it is a phase of the wave function, it could be 
related to physical effects and, in some cases, has a con- 
nection with topological numbers. The quantum Hall ef- 
fect in Bloch electron systems is described as an adiabatic 
charge transport whose process is closed in a parameter 
space. A Berry phase is generated and the quantized 
Hall conductivity is written by the Berry phase as well 
as the Chern number . Recently, it has been 

pointed out that quasiparticles in the vortex lattice of 
d x 2_ y 2-w&ve superconductors are in the Bloch states [^j. 
The spin quantum Hall effect occurs and its conductivity 
is written by a Chern number ||. Then, one can expect 
that the spin Hall conductivity in the vortex state is writ- 
ten by a Berry phase, when the effect can be described 
as an adiabatic spin transport of a closed process. 

In this paper, we discuss Bloch quasiparticles in the 
vortex state of p-wave superfluid 3 He in two dimensions 
(2D). We consider a magnetic field with a weak and ho- 
mogeneous gradient. Such a field can not be introduced 
in superconductors due to the Meissner effect. The mag- 
netic field couples to spin through the Zeeman term and 
does not to orbital currents because of the neutrality of 
the superfluid. The spin Hall current flows as an adi- 
abatic spin transport. Its conductivity is written by a 
Chern number and quantized when an excitation gap ex- 
ists. We show that the conductivity is closely related to 
the Berry phase. We also point out that the results have 
some similarity to the adiabatic pumping in Bloch elec- 
trons [flO| and the spontaneous polarization in the crys- 
talline dielectrics We set % = c = /Zb = 1, where /j,b 
is the Bohr magneton. 



II. SUPERFLUID HELIUM 3 IN TWO 
DIMENSIONS 

Let ^ a (x) stands for the Fermion field with spin a =| 
, |. The mean field Hamiltonian for fermionic superfluid 
(or superconductors) in D-dimcnsional space is written 
with the gap matrix (x, y) as 



H MF = / d D xyji( X )e(p)Mx) 



(1) 



d u xd 



A a/9 (x,y)V£(x)^J(y)+/i.c. 



P 2 



Pf 



2m 



where p = — zV and the repeated Greek indices are 
summed up. We may consider the Fourier transform of 
the gap matrix in terms of the relative coordinate of x 
and y, i.e. 



A a/3 (x,y) 



P 



!(x - y>P A a/J (r,P) 



(2) 



where r = (x + y)/2. 

Superfluid 3 He is in the spin triplet p-wave states [ p"2"| . 
In general, the gap function for spin triplet pairing is 
A(r, p) = i<jy<j ■ d(r, p), where d(r, — p) = d(r, p) is 
a three-dimensional vector in the spin space. In p-wave 
states, the magnitude of the relative angular momentum 
of the Cooper pair 1 = 1 and d vectors have a linear 
dependence on p. It is well known that the three phases 
are observed in the superfluid 3 He (A, B, and Ai phases) 
|l2| . Those phases are represented by the different d 
vectors, respectively. 

We consider a 2D system. To realize 2D, we introduce 
a strong confine potential along z-axis to avoid quasipar- 
ticle excitations along z-axis. The boundary effect intro- 
duced by the confine potential locks the relative angular 
momentum of all the Cooper pairs in the same direc- 
tion along z-axis [12] . Then, we take l z = 1 (l z : the 
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z-component of the angular momentum) in the whole re- 
gion. The direction of the d vector becomes parallel to 
the angular momentum (i.e. d//e z ) because of the ex- 
istence of the magnetic dipole interaction which couples 
spin and orbit Jj"2"| . Then, the d vector in our situation 
is 



d(r,p) = e z (f>{r)(p x +ip y ). 



(3) 



This state corresponds to the A-phasc |12| . Since the di- 
rection of the d vector and the relative angular momen- 
tum are frozen, we may neglect the textures and coreless 
vortices fL^] . The Hamiltonian for 3 He-A is 

Hmf = / d 2 x^(x)e(p)V' Q (x) 



+ / d 2 xd 2 y [A A (x,y)V{(x)Vj(y) + h.c] , (4) 



A A (x,y) = -Tr^Afoy) 



III. VORTEX STATE IN HELIUM 3 A-PHASE 
WITH A ROTATION 

It is well known that rotating superfluid is direct anal- 
ogy of type-II superconductors, and actually the vortex 
states in superfluid 3 He are detected by the experiments 
JL3| . Then, we consider superfluid in a container that 
rotates around z-axis with an angular velocity Q. Here- 
after, we use the rotating frame which is fixed on the 
container. In the rotating frame, Hmf is transformed as 



MF 



H = Hmf O • L, 



d 2 ^i( X ) 



e(p — mR) 



+ (pairing terms), 
R = Jl x x, 



(5) 



where L is the total angular momentum of Fermions 
jl4],[l5| . The kinetic energy for the quasiparticle is trans- 
formed as e(p) — > e(p — mR) — mR 2 /2. We consider 
51 ~ lrad/s and |r| < ro ~ 1mm (ro- the radius of 
the container) p3| , and we can neglect — mR 2 /2 term. 
Otherwise, we can cancel out this term by introducing 
a parabolic trap |lq] . To avoid this term is essential to 
introduce the translational invariance we will discuss be- 
low. 

Then, the one to one correspondence can be seen be- 
tween our system and charged superfluid in a magnetic 
field with an infinite London penetration depth, i.e. the 
strongly type-II superconductors. The vector field R cor- 
responds to a "vector potential A" and the Fermion mass 
m corresponds to "the electric charge e" . And "the mag- 
netic field" is 20 = V x R. 



Let us consider a vortex state, i.e. 57 > O c i and set 
up a square vortex lattice. We would like to note that 
our discussion is applicable to other types of lattices. In 
the vortex state, the gap function has the singular phase 
y(x) which satisfies, 



A A (x,y) = A A (x,y)e-*W*> + *M], 
V x W(x) = 2ne z V] S 2 (x - r<), 



(G) 



where A A (x, y) is the gauge invariant part of the gap 
function, = (e x Zj + e y n{)a with integers U and ni is 
the i-th lattice point, and e x and e y are the unit vector 
of the Cartesian coordinate in the rotating frame. When 
Q ~ lrad/s the vortex lattice constant a = y/n/mCl ~ 
10~ 2 cm. From Eq. (@) , the Hamiltonian density operator 
can be written in the Nambu representation as 

W(p,x,y) (7) 
/ e(p-mR)(S(x-y) A A (x, y ) e -*^( x )+v(y)] \ 
~ V-^A( x .y) el[v(x)+¥,(y)1 ~ e (p + mR)5(x - y) /' 

The Bogoliubov-de Gennes (BdG) equation |L7| is 

J d 2 yW(p,x, y )$ £ (y) = E* B {x), 

$ £ (x) = ([/ £ (x),-^(x)) T , (8) 

where, 

E 

WW = E [-^(x)7£T + UeWvU > 

E 

and 7I and ^Ea are the creation and annihilation op- 
erator of the Bogoliubov quasiparticles, respectively. 

Let us discuss the periodicity of the system |1S L The 
multivalued phase field y(x) which satisfies Eq.(6|) has 
an ambiguity for deformations which does not change the 
topology of its configuration, i.e. the ambiguity remains 
in terms of the gauge degrees of freedom. So, we may 
take a constraint 



ip{x + e x a) = <p(x) - ae x 
<p(x + e y a) — <^(x) — ae y 



Obviously, it is consistent with Eq. 
fine a translation operator 



mR 
mR 



T$ T = exp [iSr ■ (p + mRr 3 )] 



(9) 

Then, let us de- 
(10) 



which is the direct analogy of the magnetic translation 
operator. The symbol T3 denotes the third Pauli matrix 
in the Nambu (particle-hole) space. The coordinates are 
translated by T$ r as x — * x + Sr and y — > y + Sr. Then, 
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one can see easily that the operator T exa and T Gya com- 
mute with W(p,x, y), but does not commute with each 
other. 

We define a unit cell, in which a "unit flux" 2n/m 
penetrates. It is a direct analogy of the magnetic unit 
cell, where the magnetic unit flux 27r/e penetrates. A 
vortex has a "half unit flux" ir/m and two vortices are 
contained in a unit cell. Assume that there are even 
numbers of vortices, and one may choose the unit cell as 
Fig. 1. Consider the translations in terms of the cell, 



T en .a+e y a and T e 



L e x a—e y aj 

One can see easily that the operators satisfy 



where d = v2a. 



[W(p,x,y),T 5r 



0. 



(11) 



Therefore, the eigenstate o/7i(p,x, y) are in the Bloch 
state, i.e. 



*k(x) 



C "k(x), 



(12) 



where k is in the Brillouin zone (BZ), —iz/d < (k x , k y ) < 
ir/d. Here, we omit the band index. Define 



H k (x,y)^ e - ikx H(p,x, y ) e ik y. 
From Eq.(@) 

Wk(x,y) -H(p + k,x,y) e - i;k -( x -y) 



(13) 



(14) 



Then, from Eq.f| and @, one can see that the function 
«k(x) satisfies 



J d 2 yH k (x,y)u k (y) = E k u k (x), 



(15) 



and its translation in terms of the unit cell satisfies a 
generalized Bloch condition H 



Uk(x + e' x d) = exp [ide x • rnRr^} itk(x), 
itk(x + e' y d) = cxp [ide' y ■ mRr 3 ] «k(x). 



(16) 



del 



del 



ae 




ae 



FIG. 1. The unit cells surrounded by green dotted 
lines. Black dots denote the vortices. 



The spectrum for lattice quasiparticles in the square 
vortex array of the p x + ip y -w&ve superconductors has 
been investigated by using the singular gauge transfor- 
mation p[Jl9|]. A zero energy state and gaps around it 
are found. In our situation, the continuum system is 
considered. There might be some problems to apply the 
singular gauge transformation approach to the contin- 
uum systems. It was pointed out that the quasiparticle 
spectrum depends on the choice of the singular gauge 
transformations |2l]]. Here, we assume the existence of 
an excitation gap. 



IV. THE ADIABATIC PROCESS AND THE 
BERRY PHASE 



Let us introduce a magnetic field, which is directed to 
z-axis and has a homogeneous gradient in the rotating 
frame, which is written B z (x) — x • VB Z and V£> 2 is 
a constant vector. The field will be a driving force of 
the spin transport. For a moment, we consider in the 
Lagrange formalism. In supcrfluid, the magnetic field 
couples to spin through Zeeman term and does not cou- 
ple to orbital currents. Then, the Lagrangian is written 
in the form 

C = J d 2 x^(x) {id/do - (x • VB z /2)} *(x) 

- J d 2 xd 2 y^^)H(p^,y)^(y), 
*(x)^ T (x),^(x)) T (17) 

We consider a phase transformation of Eq.(|l7j), 

*(x) -> exp[-itx • VB»/2]$(x). (18) 



Then, the term — (x- VB z /2) is absorbed and the Hamil- 
tonian density operator is transformed as 



where, 



H(p, x, y) -> H(f> - f (*), x, y), 



f(i) = WB z /2. 



(19) 



(20) 



By using the analogy of U(l) (electromagnetic) gauge 
theory, we may regard f(t) as a vector potential that 
couples to the spin current, since it is introduced by the 
local spin rotation Eq. fll8| ) . We assume that f (i) changes 
adiabatically, i.e., |V-B Z | << 1. For simplicity, we write 



«(p-f(t),x,y) = «(f,x,y). 



(21) 



Then, we solve a time-dependent equation of motion with 
the adiabatic parameter f(i) 
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d f 
i— *(t,x) = J d 2 yU{t,^y)^(t,y) 



(22) 



We use the adiabatic approximation and an eigenvalue 
equation at fixed t is 



J d 2 yH{t, X ,y)<P E{t) {t,y) = E(t)$ B(t) (t,x) 



(23) 



Obviously, it is equivalent to the BdG equation (g) at 
t = 0. The Hamiltonian 7i(t, x, y) has a spatial periodic- 
ity as well as Ti(t = 0, x, y) because VB Z is homogeneous. 
Then, eigensolutions are written in the Bloch form, i.e., 
$k(*,x) = e 4kx w k (i,x) (See, Eq. ©). The function 
Mk(i,x) obeys the equation, 



J d 2 yHk{t 7 x, y)u k (t, y) = E k (t)u k {t, x), (24) 
W k (t,x,y)=W k _ f(t) (x,y), (25) 

and hence, 

u k (t,x) = u k _ f(t) (x). (26) 
The solution of Eq. (B2fl in the adiabatic approximation 



*k(*,x), 



*k(i,x) = exp 

Tt(*)=*jf df/*k(f 



i I dt'(E k (t')+ lk (t')) 
o 



<9f' 



(27) 



The reciprocal lattice vector for the square vortex array 
is written G = (le' x + ne' y )(2-K / d) with integers I and n. 
We note that it is possible to compactify the Hamilto- 
nian as W k (i,x, y) ~ Hk+G(t, x, y), because they give 
the equivalent eigensolutions. Also we note that the pa- 
rameter f (t) varies on the Brillouin zone, (See, Eq. (p5|)). 
Therefore, when f(t)//G we have a period T for a closed 
loop on the parameter space pCfl. For example, 



T = A7r/(\VB z \d), 



(28) 



for f(i)//e' x ,e' y and the Berry phases are defined as 

Jq dtjk(t) for each cases Q. We introduce the Berry 
connection, a(k) = (u k |V k |u k ), which is a gauge field 
defined on the parameter space Q. By using Eq. (26), 
the Berry phases for the parameter f (t)/ /e' x , e' y are writ- 
ten as 



and 



dk x a x (k), 



dk y a v (k), 



(29) 



(30) 



respectively. 

When f(t)//G, we could write down the Berry phase 

as 



T f = if dk-a(k), 

JC(f) 



(31) 



where C(f) is a closed loop on which f(£) moves. In gen- 
eral, the Berry phase depends on C(f) [0. 



V. SPIN QUANTUM HALL EFFECT AND THE 
BERRY PHASE 



Let us calculate the spin current. In 3 He-A, the sys- 
tem is invariant under the spin rotation around z-axis 
$ k (t) -> e ,e $ k (t) and H(f,x,y) -> e i9 H{t, x, y)e~ ie . 
The spin current j s is defined by the spin conservation 
law p2[ , i.e., p s + V • j s = 0, where p s is the spin density 
_(2tt)V(x) = (l/2) E„<o/ B z d2fcvI/ ikW*nk(x) and we 
introduce the band index n. The label denotes the zero 
energy. As we mentioned before, we assume an excita- 
tion gap, i.e. there are no partially filled bands. Then, 
the response of the spin current for the uniform field f (t) 
is 



<n*)} 



It. 



n<0 



E 



d 2 k 

LiZ W) 

d 2 k 

LiZ (27T) 2 



■2 

Unk(0 



1 



n<0 

-a*, [VB, x e z ] , 



[r,W(t)] 
5u nk (t) 



*nk(*) 



/l.C. 



(«) 

Ch ' 



where 



N. 



n<0 



(n) 



Ch 



d 2 k 



BZ 



2t™ 



t [V k x a„(k) 



(32) 



(33) 



is the Chern number for the n-th band and a„ (k) is equiv- 
alent to the Berry connection for the n-th band. The 
detailed calculation is written, for example, in Ref. JtJ. 

The Chern number takes integer. The reason is based 
on the fact that a n (k) is defined on the torus (the BZ) 
and the Chern number becomes finite if and only if a k has 
a non-trivial topology. The nature of the Chern number 
has been discussed in detail in Ref. ||. 

Then, Eq. d32| ) shows that a stationary spin Hall cur- 
rent flows as an adiabatic spin transport and its conduc- 
tivity is quantized as an integer multiple of 1/8tt. The 
same result for the conductivity in the vortex state of 
d-wave superconductors have been obtained Q . The dis- 
crete conductance change is expected to occur when one 
varies fl and pp . 
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Here, we calculated the expectation value of the total 
spin current directly by using the adiabatic approxima- 
tion and obtained the spin Hall conductivity. One can 
see the fact in the calculations of Eq. (32) that this ap- 



proach is equavalcnt to calculate the Kubo formular for 
the spin Hall conductivity argued in Ref. ||. 

Finally, we show that a s could be written in terms of 



the Berry phase [[l]|5j,g,|7|] ■ By using the Stokes' theorem, 
we have the relation (See, Eq.@, Eq. (§(]), Eq. (§|) and 
Eq. (H)), 



-xy 



—y 

16tt 2 ^ 



(34) 



n<0 



dk 



dT n v (k x ) 



dk 



VI. RELATIONS TO THE OTHER ARGUMENTS 

The effect have some similarity to the adiabatic pump- 
ing which is originally argued by Thouless and discussed 
actively at present |10| . In pumping, an adiabatic ac 
perturbation yields a dc current, and the charge trans- 
fer per a cycle is independent of the period of the per- 
turbation. The charge transfer is quantized when the 
ac perturbation is commensurate with the lattice in ID. 
As we mentioned before, the Hamiltonian in our system 
Hk(t) is compactified and moved periodically by the adi- 
abatic parameter t(t). Then, the Hamiltonian changes 
ac-like, and the change yields a dc spin Hall current. To 
make a correspondence to the Thouless' arguments, one 
calculates a spin transfer per period T. Assume that 
i(t) = e' y t\V B z \/2. The spin Hall current flows along x'- 
axis (See, Fig. 1) and the spin transfer per the boundary 
of the unit cell along y'-axis is 



r T 

AS z = d dt(f x ,(t))=-Y / - 

J ° n<0 



iV Ch 



(35) 



(see, Eq. (|32J) and (|3J)). The result does not de- 

pend on T . It comes from the fact that both of the 
magnitude of the quantized current and T^ 1 are propor- 
tional to |V.B Z |. We emphasize that the spin transfer 
is quantized, i.e., the integral spin transfer occurs. The 
result is analogous to the Thouless' result M]. 

The value AS z /d corresponds to the magnetization 
change per the period. From Eq. (34), the magneti- 
zation change is written by the Berry phase. Then, the 
present result is also similar to the spontaneous polar- 
ization of crystalline dielectrics , which is written by the 
Berry phase introduced by a closed adiabatic change of 
the Kohn-Sham potential jnj . 

Essentially, the similarity comes from the fact that the 
effects argued here are caused by the closed adiabatic 



change in the Bloch states with the finite energy gap. A 
parallel discussion for the present arguments have been 
made in the Bloch electron systems in the presence of the 
electromagnetic field with respect to the charge transport 

0- 



VII. SUMMARY AND DISCUSSIONS 

In summary, we consider Bloch quasiparticles in a vor- 
tex state of superfluid 3 He-A in 2D with a rotation along 
z-axis. A magnetic field is along z-axis with a weak 
homogeneous gradient in the rotating frame. The field 
could be represented by an adiabatically changing vector 
potential which couples to the spin current. The adia- 
batic process is defined on a closed loop in the parameter 
space (the Brillouin Zone) and generates a Berry phase. 
The spin Hall current flows in the process. We calcu- 
lated the expectation value of the total spin current di- 
rectly by using the adiabatic approximation and obtain 
the spin Hall conductivity. This approach is equavalent 
to calculate the Kubo formular for the spin Hall conduc- 
tivity The conductivity is represented by the Chern 
number and quantized when the quasiparticle has an ex- 
citation gap as that in the cJ-wave vortex state ||. We 
have shown that the spin Hall conductivity is written 
by the Berry phase. The spin transfer per a cycle per 
the boundary of the unit cell is quantized and related 
to the Berry phase. The results remind us the adiabatic 
pumping, which is introduced by Thouless with respect 
to the charge transport . The result is also similar to 
the relation between the spontaneous polarization and 
the Berry phase in the crystalline dielectrics |ll]] . Essen- 
tially, the similarity comes from the fact that the effects 
argued here are caused by the closed adiabatic change in 
the Bloch states with the finite energy gap. With respect 
to the charge transport, a parallel discussion have been 
made in the Bloch electron systems in the presence of the 
electromagnetic field . 

As we mentioned before, the spin quantum Hall effect 
in the vortex state of c4;2_ y 2-wave superconductor has 
been pointed out |J, but it seems to have some diffi- 
culty to make a parallel discussion in the superconduc- 
tors. Because of the Meissner effect, it is not possible to 
have a magnetic field with a finite homogeneous gradient 
which is essential to define the adiabatic process on the 
closed loop in the parameter space. The vortex states 
in 3 He-A are suitable for our arguments because 3 He is 
the fermionic superfluid in which the spin current is well 
defined, i.e. the spin rotation symmetry around z-axis is 
remained. In contrast to the c? a; 2„ l) 2-wave state, the spin 
quantum Hall effect occurs spontaneously in 3 He-A , i.e. 
one obtains a quantized spin Hall conductivity to cal- 
culate the Kubo formular in the absence of the vortices. 
The effect comes from the broken time reversal symmetry 
and the broken parity in the orbital part of the pairing 
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symmetry p2[ . But, the system does not have the fi- 
nite spatial periodicity and we can not make a parallel 
discussion also in this case. 

We would like to comment on the fact that, in our ar- 
gument, the orbital part of the pairing symmetry is not 
crucial as long as the quasiparticle spectrum in the vor- 
tex state has an excitation gap. The spin part is crucial 
because the spin rotational symmetry is needed to obtain 
well defined spin currents. 

Several authors have made their efforts to find out a 
way to measure spin transport p3[ |. Some experimental 
techniques to detect spin transfer is highly desirable. 

The authors are grateful to K. Maki, M. Sato, Z. 
Tesanovic and F. Zhou for useful discussions. 
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